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Atom interferometers covering macroscopic domains of space-time are a spectacular man- 
ifestation of the wave nature of matter. Due to their unique coherence properties, Bose- 
Einstein condensates are ideal sources for an atom interferometer in extended free fall. In 
this paper we report on the realization of an asymmetric Mach-Zehndcr interferometer op- 
erated with a Bose-Einstein condensate in microgravity. The resulting interference pattern 
is similar to the one in the far-field of a double-slit and shows a linear scaling with the time 
the wave packets expand. We employ delta-kick cooling in order to enhance the signal and 
extend our atom interferometer. Our experiments demonstrate the high potential of inter- 
ferometers operated with quantum gases for probing the fundamental concepts of quantum 
mechanics and general relativity. 
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Quantum theory [JJ and general relativity (2; are two pillars of modern physics and successfully 
describe phenomena of the micro- and the macro-cosmos respectively. So far they have resisted 
any attempt of complete unification and quantum gravity [3] is generally considered the Holy Grail 
of physics. Experimental tests of gravity [I] with matter waves [5] started as early as 1975 with 
neutrons (HE]. Today, atom interferometers (AI) [8] offer new opportunities to probe the interface 
of these fundamentally disparate descriptions of nature. The coherent evolution of quantum objects 
delocalized in space-time [S], the verification of the Einstein principle of equivalence with quantum 
objects [10] and the detection of gravitational waves (TTJ constitute only three of many timely 
quests motivating experiments with AI in extended free fall. The overarching aim is to enhance 
the sensitivity of these devices, which increases linearly with the momentum difference between the 
two matter waves |12] emerging from a beam splitter and quadratically with the time of free fall 
as experienced in fountains [TQl [13] , drop towers p3] , parabolic flights [15] and space [16] . These 
scaling laws imply constraints with respect to the atomic source. Thanks to their slow spreading and 
their excellent mode properties, Bose-Einstein condensates (BECs) \17\ [TH] represent a promising 
source [§] for high-resolution interferometers |19H21| . 

In this letter, we report on the demonstration of a BEC-interferometer in microgravity. Taking 
advantage of the extended free fall provided at the drop-tower of the Center of Applied Space 
Technology and Microgravity (ZARM) in Bremen, we have been able to coherently split a BEC 
consisting of about 10 87 Rb atoms and to separate the emerging wave packets over macroscopic 
scales in time and space. The interferometer extends over more than half a second and covers 
distances of millimeters, exceeding the width of the condensate by an order of magnitude. By 
applying delta-kick cooling (DKC) [22-23] we have been able to reduce the expansion and to 
enhance the signal at longer interferometry times. 

Our experiments are performed with an asymmetric Mach-Zehnder interferometer (AMZI) pJU 
[25] shown in Fig. 1. Here we display the temporal evolution of the atomic density distribution of the 
BEC-interferometer for an experiment on ground [Fig. 1(a)] and the corresponding experimental 
sequence for forming the interferometer [Fig. 1(b)]. A macroscopic wave packet is coherently split, 
redirected and brought to a partial overlap by successive Bragg scattering at moving light-crystals 
[26H28] , They are generated by pulses of two counter-propagating laser beams separated by the 
two-photon recoil energy of 15 kHz and detuned from the F=2 — > F=3 transition of the D2 line 
of 87 Rb by 800 MHz to reduce spontaneous scattering. There exists a close analogy to the Young 
double-slit experiment [Fig. 1(c)], where one pair of overlapping BECs plays the role of a pair of 
coherent light waves emanating from two slits separated by a distance d. Similar to the resulting 
interference pattern in the far-field of the double-slit, the fringe spacing in our expanding cloud, 
being the distance between two local maxima of the density, increases with the total expansion 
time T ex of the BEC and is inversely proportional to the displacement d of the two clouds. 

Figure 2 illustrates our experiments in microgravity performed at the drop tower. We show 
the complete temporal sequence [Fig. 2(a)] which differs from the previous experiments with the 
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Figure 1. Temporal AMZI for a BEC based on Bragg scattering at a light grating: experimental images on 
ground (a), schematic sequence (b) and analogy to the Young double-slit experiment (c). The evolution of 
the BEC and the AMZI is visualized by a series of absorption images (a) of the atomic densities separated 
by 1 ms. The incomplete transfer is a consequence of the larger mean field energy of the BEC necessary 
for the ground experiment. The interferometer starts at the time T after the release of the BEC, when a 
7r/2-pulse (b) made out of two counter-propagating light beams of frequency u and lo + S creates a coherent 
superposition of two wave packets which drift apart with the two-photon recoil velocity v rec = 11.8 mm/s. 
After T they are re-directed by a 7r-pulse and partially recombined after T — ST by a second 7r/2-pulse. A 
nonzero value of ST leads to a spatial interference pattern, which we record after r = 53 ms in free fall. 
Similar to the far-field pattern observed in the Young double-slit experiment (c), the fringe spacing A scales 
linearly with the time of expansion T ex = Tq + 2T — ST + r and is inversely proportional to the separation 
d = v rec ST of the wave packets. The scaling factor is the ratio of Planck's constant h and the mass mub of 
the Rubidium atoms. 



apparatus [2] in three important features: (i) We employ DKC to reduce the expansion during 
the free fall by briefly (2 ms) switching on the trap with frequencies of (10, 22, 27) Hz generated 
by the atom chip 30 ms after the release, (ii) In order to eliminate detrimental effects of residual 
magnetic fields we transfer the BEC into the non- magnetic state \F = 2, tuf = 0) by coupling the 
Zeeman levels with a chirped radio- frequency pulse (adiabatic rapid passage), (iii) At the time 
To after the release of the BEC, we implement the sequence of the AMZI outlined in Fig. 1 and 
detect the interference pattern at T ex after the release using absorption imaging with a single laser 
pulse as illustrated in Fig. 2(b). In Fig. 2(c) we show typical images of the interfering BECs and 
the corresponding column density profiles for two different values of T ex . 
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Figure 2. Mach-Zehnder interferometry of a BEC in microgravity as realized in the ZARM drop-tower in 
Bremen (a) where absorption imaging (b) brings out the interference fringes (c). The preparatory experimen- 
tal sequence (a) includes capturing cold atoms in a magneto-optical trap (MOT), loading a Ioffe-Pritchard 
trap, creating a BEC and applying the DKC followed by the adiabatic rapid passage (ARP). The remaining 
time before the capture of the capsule at the bottom of the tower is used for atom interferometry (AI) and 
imaging of the atoms. The AMZI below the atom chip (top plane of b) is formed by scattering the BEC 
off moving Bragg gratings generated by two counter-propagating laser beams (red arrows directed along 
the y-axis), resulting in two pairs of interfering BECs. A resonant laser beam propagating along the x-axis 
projects the shadow of the BEC onto a CCD camera. Typical interference patterns and the corresponding 
column densities (c) are shown for T ex of 180 ms and 260 ms with corresponding fringe spacing of 75 pm 
and 107 pm. 

Figure 3 summarizes the central results of our paper on probing the coherent evolution of a 
BEC with an interferometer in extended free fall. It shows the spatial period of the observed 
fringe pattern [Fig. 3(a)] as a function of the expansion time T ex , and the contrast [Fig. 3(b)] 
observed at the exit ports of the AI for increasing values 2T — ST of the time the BEC spends in 
the interferometer. Moreover, in Fig. 3(a) we confront the experimental results (blue circles, red 
squares, black triangles) with the corresponding theoretical predictions (solid blue and red lines). 
The solid blue line originates from a model based on the scaling approach [29] and describes the 
interference pattern of two condensates initially separated by a distance d, which start to expand 
and eventually overlap. Their initial shape is derived from a detailed numerical model of our 
magnetic chip trap. For large time scales the observed fringe spacing (blue dots) shows a linear 
increase with T ex in full accordance with our model and with the linear far-field prediction (dash- 
dotted blue curve) of the double-slit. We emphasize, that our microgravity experiments operate 
deep in the linear regime and the non-linear behavior typical for the near-field combined with the 
non- linear evolution of the BEC occurs only at very short times (<30 ms). The linear scaling of 
the fringe pattern confirms the unperturbed evolution of the BEC during extended free fall. 
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Figure 3. Fringe spacing (a) of two interfering BECs observed at each exit port of an AMZI with (red 
squares, black triangles, solid red line) and without (blue dots, solid blue line) delta-kick cooling (DKC) as 
a function of the expansion time T ex , and contrast (b) in its dependence on the time 2T — ST spent in the 
interferometer. Our data (blue dots, red boxes, black triangles) agrees well with our models of free expansion 
(solid blue line) and DKC (solid red line). In the case of DKC, the temporal asymmetry ST was increased 
from 1 ms to 2.5 ms in order to preserve the number of visible fringes across the smaller BEC by rescaling 
the fringes, which decreased the slope in (a). With and without DKC, the observed linear dependence has 
the slope of the approximation of the double-slit for point particles (dash-dotted lines). An offset arises 
from the initial size and the nonlinear expansion of the BEC due to mean-field or from the focusing in DKC 
which shifts the apparent location of the double-slit to earlier or later times, respectively. The high contrast 
of the interference fringes (b), obtained by absorption imaging, typically exceeds 40% but fades away for 
increasing T. Error bars depict 1-tr-confidence bounds of the fitted parameters. 



The expansion rate of the BEC due to the mean-field energy is a limiting factor for extending the 
interferometer to even longer time scales. It can be reduced by DKC, which acts on the BEC like a 
three-dimensional lens. Indeed, in our experiments DKC realized with the atom chip eliminates a 
substantial part of the kinetic energy of the BEC, giving rise to an effective temperature of about 
InK. The method allows us to extend the observation of the free evolution of the BEC and was 
tested with our AMZI. The experimental observations (red squares) of the fringe spacing agree well 
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with the theoretical predictions (solid red line) for a double-slit experiment with delta-kick cooled 
atoms. In order to reach even longer times (black triangles), we have adjusted the detection time 
r such that the patterns of the two exit ports overlap, thus increasing the absorption signal. 

As shown in Fig. 3(b), we observe a contrast of more than 40%, even at times 2T — 5T as 
large as half a second. However, then the contrast decreases with the time over which the wave 
packets are separated, and generally with the expansion time of the BEC. The observed reduction 
is non-exponential in time and uniform over the cloud. In this respect the asymmetric interfer- 
ometer puts more severe constraints on the set up than the symmetric one. However, it allows 
us to analyze various effects perturbing the interferometer. A preliminary analysis shows that the 
reduction may be due to an imperfect alignment of the beam splitters, inhomogeneous wavefronts 
or disturbances resulting from a slight capsule rotation. A more detailed discussion is subject to 
future investigations. 

In conclusion, we have followed the evolution of the temporal coherence of an asymmetric 
Mach-Zehnder interferometer with a BEC in microgravity, which is analogous to Young's double- 
slit experiment with a gigantic matter wave packet being in free fall for nearly a second. All 
preparatory steps necessary for high-precision interferometry are implemented with the help of 
a robust chip-based BEC source. Our device represents a unique testbed for exploring atom 
interferometry with novel states of matter in extended free fall. In particular, it allows us to 
test tools of atom optics, such as the precise mode control of BEC with DKC at energy scales 
approaching pK temperatures. These concepts are essential for high-resolution measurements 
both in fountains and in microgravity. Moreover, the free fall of our apparatus results in a rather 
compact interferometer. Indeed, a fountain-like interferometer would need more than hundred 
meters to achieve the time of free fall provided by the catapult at the ZARM. Last but not least, 
our experiment paves the way for quantum tests of the weak equivalence principle on sounding 
rockets, the International Space Station and satellite missions [30] such as STE-QUEST. 
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Centre for Quantum Engineering and Space-Time Research. 



The Authors contributed equally 



1 To whom correspondence should be addressed: rasel@iqo.uni-hannover.de 



[1] D. Bohm, Quantum Theory (Prentice Hall, Englewood Cliffs, 1951). 

[2] C. W. Misner, K. S. Thorne, J. A. Wheeler, J. Wheeler, and K. Thorne, Gravitation, first edition ed. 

(W. H. Freeman, San Francisco, 1973). 
[3] C. Kiefer, Quantum Gravity, 2nd ed. (Oxford University Press,, Oxford, NY, 2007). 



7 



[4] 
[5] 
[6] 
[7] 

[8] 
[9] 
[10] 



C. M. Will, |Living Rev. Relativity 9 (20 06) 

P. R. Berman, Atom Inter ferometry (Academic Press, San Diego, CA, 1997). 



R. Colella, A. W. Overhauser, and S. A. Werner, Physical Review Letters 34, 1472 (1975) 



H. Rauch and S. Werner, Neutron Interferometry: Lessons in Experimental Quantum Mechanics (Ox- 
ford University Press, Oxford, NY, 2000). 



A. D. Cronin, J. Schmicdmayer, and D. E. Pritchard, Reviews of Modern Physics 81, 1051 (2009) 



B. Lamine, R. Herve, A. Lambrecht, and S. Reynaud, |Physical Review Letters 96, 050405 (2006) 



S. Dimopoulos, P. W. Graham, J. M. Hogan, and M. A. Kasevich, Physical Review Letters 98, 111102 
|(2007)[ 



S. Dimopoulos, P. W. Graham, J. M. Hogan, M. A. Kasevich, and S. Rajendran, Physical Review D 
787122002 (2008)] 



[11] 
[12] 

[13] H. Muller, S.-w. Chiow, S. Herrmann, S. Chu, and K.-Y. Chung, |Physical Review Letters 100, 031101 
[14] 



S.-w. Chiow, T. Kovachy, H.-C. Chien, and M. A. Kasevich, physical Review Letters 107 , 130403 
|(2011)[ 



(2008) 



[15] 

[16] 
[17] 
[18] 
[19] 

[20] 



T. van Zoest, N. Gaaloul, Y. Singh, H. Ahlers, W. Herr, S. T. Seidel, W. Ertmer, E. Rasel, 
M. Eckart, E. Kajari, S. Arnold, G. Nandi, W. P. Schleich, R. Walser, A. Vogel, K. Sengstock, 
K. Bongs, W. Lewoczko-Adamczyk, M. Schiemangk, T. Schuldt, A. Peters, T. Konemann, H. Miintinga, 
C. Lammerzahl, H. Dittus, T. Steinmetz, T. W. Hansch, and J. Reichel, [Science 328, 1540 (2010)| 
R. Geiger, V. Menoret, G. Stern, N. Zahzam, P. Cheinet, B. Battelier, A. Villing, F. Moron, M. Lours, 
Y. Bidel, A. Bresson, A. Landragin, and P. Bouyer, |Nat Commun 2, 474 (201 1)[ 
E. Arimondo, W. Ertmer, W. Schleich, and E. M. Rasel (IOS Press, Bologna, Italy, 2009). 



E. A. Cornell and C. E. Wieman, [Reviews of Modern Physics 74, 875 ( 2002) 



W. Ketterle, |Reviews of Modern Physics 74, 113l~(2 002) 



Y. Torii, Y. Suzuki, M. Kozuma, T. Sugiura, T. Kuga, L. Deng, and E. W. Hagley, Physical Review 
|A 61, 041602 (2000)| 

J. E. Simsarian, J. Denschlag, M. Edwards, C. W. Clark, L. Deng, E. W. Hagley, K. Helmerson, S. L. 
Rolston, and W. D. Phillips, [Physical Review Letters 85, 2040 (2000) 



[21] J. E. Debs, P. A. Altin, T. H. Barter, D. Doring, G. R. Dennis, G. McDonald, R. P. Anderson, J. D. 



Close, and N. P. Robins, Physical Review A 84, 033610 (2011) 



[22] 
[23] 
[24] 
[25] 
[26] 

[27] 

[28] 
[29] 
[30] 
[31] 



S. Chu, J. E. Bjorkholm, A. Ashkin, J. P. Gordon, and L. W. Hollberg, [Optics Letters 11, 73 (1 986) 



H. Ammann and N. Christensen, Physical Review Letters 78, 2088 (1997) 



M. Morinaga, I. Bouchoule, J.-C. Karam, and C. Salomon, Physical Review Letters 83, 4037 (1999) 



H. Rauch, H. Wolwitsch, H. Kaiser, R. Clothier, and S. A. Werner, Physical Review A 53, 902 (1996) 



M. Kozuma, L. Deng, E. W. Hagley, J. Wen, R. Lutwak, K. Helmerson, S. L. Rolston, and W. D. 
Phillips, physical Review Letters 82, 871 (1999)[ 

Y.-J. Wang, D. Z. Anderson, V. M. Bright, E. A. Cornell, Q. Diot, T. Kishimoto, M. Prentiss, R. A. 
Saravanan, S. R. Segal, and S. Wu, [Physical Review Letters 94, 090405 (2005)) 



K. Nakagawa and M. Horikoshi, Journal of Physics: Conference Series 185, 012031 (2009) 
See Appendix for expansion theory and image analysis. 

ESA, The space missions selected by ESA's Cosmic Vision 2020 22 are listed at http:/ /sci.esa.int. 
J. C. Camparo and R. P. Frueholz, Journal of Physics B: Atomic and Molecular Physics 17, 4169 
(1984) 



s 



[32] Y. Castin and R. Dum, Physical Review Letters 77, 5315 (1996) 



[33] P. Storey and M. Olshanii, Physical Review A 62, 033604 (2000) 



[34] G. Nandi, R. Walser, E. Kajari, and W. P. Schleich, Physical Review A 76, 063617 (2007) 



Al 



Appendix 

Supplementary Material for 

Interferometry with Bose-Einstein Condensates 

in Microgravity 

MATERIALS AND METHODS 

In this supplementary material, more details about key experimental steps and the theoretical 
description are given. In the experimental part we provide details on the drop tower, the generation 
of the light grating for splitting and recombining the atomic ensembles, the state preparation & 
detection and the quantitative analysis of the absorption images of fringe patterns. 

The theoretical description of the experimental results starts with a section briefly reviewing 
the scaling approach as a tool for describing the expansion dynamics of a Bose-Einstein condensate 
(BEC) for a time-dependent trapping potential. This is then applied to the asymmetric Mach- 
Zehnder interferometer (AMZI) of our experiment and all the relevant features of the interference 
fringes at the exit ports are discussed. The complementarity of the interference patterns of the two 
exit ports as well as the case of the detection of overlapping ports are discussed in two separate 
sections. 

Drop tower 

Our free-fall experiments are carried out in the drop tower of the Center of Applied Space 
Technology and Microgravity in Bremen, Germany. Inside the 146 m tall tower, a 122 m high steel 
tube 3.5 m in diameter can be evacuated to around 10 Pa within 90 minutes. After a free fall of 
4.7 s with residual accelerations in the parts per million range of the local gravity g, the sealed 
experimental capsule is caught in a 8 m high pool of polysterene balls. During the impact, the 
experiment has to sustain a deceleration of up to 50 g. A catapult mechanism housed beneath the 
tower allows us to launch a capsule upwards and thus double the free-fall time, which we plan to 
use in future experiments. 

Light-induced Bragg scattering 

Two light fields employed for Bragg scattering are generated by a single distributed feedback 
laser diode, which is stabilized with a detuning of 800 MHz to the F = 2 — > F' = 3 transition 
in the D2-line of 87 Rb. The two beams pass an acousto-optic modulator (AOM) to generate the 
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required relative offset frequency of 15 kHz. This technique provides us with precise control of the 
timing and the shape of the pulses. The light is then coupled into two polarization-maintaining 
fibers and guided to the vacuum chamber, where the two counterpropagating collimated beams 
of parallel linear polarization and diameter of 0.65 cm (FWHM) are formed. They are aligned 
parallel to the chip surface and perpendicular to the detection axis so that the interference fringes 
can be detected by means of absorption imaging. In the data presented, we use box-shaped pulses 
and Gaussian-shaped pulses. The experiments showed that the pulse shape does not influence the 
contrast nor the signal of the interference pattern in the presented data. For this reason, we did 
not distinguish between the pulse shapes in Fig. 3. 



For the preparation of the BEC in the non-magnetic state \F = 2, mp = 0) we an adiabatic rapid 
passage 5 -* (ARP). For the state transfer we apply a homogeneous static magnetic field of 11 G in 
x-direction and combine this with a 3.8 ms long RF-sweep from 7.714 MHz to 7.754 MHz. After 
the interferometer sequence a quadrupole field in Stern-Gerlach configuration is applied to split the 
atoms according to their hyperfine-state and to ensure that only atoms in the non-magnetic state 
contribute to the interference pattern. In this way we are also able to verify transfer efficiencies to 
the non- magnetic state of about 90%. 



The atomic clouds are probed with a laser beam on the F = 2 — > F' = 3 transition where the 
intensity of the detection beam is typically about 5% of the saturation intensity of the transition. 
Two images are taken in each experiment with a charge-coupled-device (CCD) camera (Hamamatsu 
C8484-15G): First the shadow of the cloud in the detection beam is recorded providing us with the 
intensity I a toms- This beam heats the atoms and expels them from the imaging region. Then the 
detection beam is imaged in the absence of any atoms yielding Ibeam- Both intensities I atoms and 
heam are corrected for the camera dark image Idark and we obtain the optical density 



State preparation 



Detection and analysis of absorption images 




(1) 



and the atomic column density 



D 



(2) 



n = — 
a 



with the resonant cross section a. 

To extract the contrast and the fringe spacing from the images of the interferometer output 
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ports, we sum the column densities along z and fit the function 



n\D = n max iD 



1 + Csin \^{y -yi) + y 



V A 



exp 



(y - yif 



T Tlmax ,1D 



2tt 



1 + Csin ( — (y - y 2 ) + (f + tt 



exp 



(y - V2f 



+ n ,iD 



(3) 



to the ID profiles with a nonlinear least squares method. Here, y\ and y<i are the centers of 
the two output ports, a y the width of their Gaussian envelopes, A the fringe spacing and C the 
contrast shown in Fig. 3(a) and Fig. 3(b) respectively, while riQ i£, accounts for an offset due to 
intensity fluctuations between the images. The phase cp would be read out for an interferometric 
measurement of e.g. inertial forces. For long observation times a tilt of the interference fringes with 
respect to the z-axis is observed, which is compensated by aligning the image before we calculate 
the column density along z. Error bars in Fig. 3 are standard deviations calculated via the fitted 
model and multiple noise realizations. 

Figure |A4| displays the decaying signal to noise ratio 



SNR = nmax ' 1D (4) 

"rms 

of the integrated atomic density with and without delta kick cooling (DKC) as a function of 
time, with u rms being the RMS of the fit residuals. Due to the spreading of the atomic wave packet, 
in the linear regime the ID atomic density decreases according to 

N „ N 

n m ax,lD = -7= f— — (5) 

vo7r<7 y V87T (aT ex + p) 

with the total atom number N, expansion rate a and an offset j3. Thus, in the case of DKC 
the SNR falls off with a slower rate and is doubled by overlapping the output ports, allowing us to 
detect the atoms after longer times T ex than in the free case. 

For longer interferometer times IT — 5T the fit routine does not converge to the expected fringe 
period, which is kept as a free parameter in the data evaluation. These data (open red squares and 
open black triangles in Fig. |A4| ) were not included in Fig. 3. Fitting the discarded data with the 
fringe period predicted by Eq. ^ yields a finite contrast. These values were confronted with the 
contrast obtained by applying the same fit routine to synthetic test signals, which were generated 
from Eq. [3] for a vanishing contrast (C = 0) with multiple noise realizations in compliance with 



the observed SNR. As depicted in the inset of Fig. A4, the contrast resulting from the analysis 
of the test data (dashed line) was indeed equal or higher than the one obtained for the discarded 
experimental data, thus confirming the loss of contrast. The dashed curve indicates that the 
minimum detectable contrast is limited by the SNR. 
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Figure A4. Signal to noise ratio (SNR) of the integrated atomic density forming the interference patterns 
as a function of the total expansion time T ex with (red squares) and without (blue dots) delta kick cooling 
(DKC) as well as for overlapping interferometer ports (black triangles). With DKC the SNR decreases 
with a slower rate than without and is further increased when both interferometer ports are detected before 
they are fully separated spatially. The scatter is caused by atom number fluctuations and by the detection 
intensity varying between individual drops due to the atoms being located in another position within the 
Gaussian profile of the laser beam. To facilitate the comparison of the three cases, the curves (blue, red, 
black) show the decay according to Eq.Q and ^ using the measured expansion rates, mean atom number 
and noise. Towards longer expansion times (open black triangles and red squares), the period of the fringe 
pattern could not be obtained by fitting the data with the period as a free parameter, so these data points 
were discarded in Fig. 3. The inset presents the observed contrast in relation to the SNR for all datapoints. 
The contrast was obtained by fitting the function n\r> given by. Eq. ([3| to the experimental data (blue 
circles, red squares and black triangles) and to synthetic data (black dashed line) generated from Eq. ^ for 
a vanishing contrast (C = 0) and noise in compliance with the observed SNR. For the synthetic data as well 
as for the discarded experimental data points the contrast was derived from fits with the fringe period fixed 
to the theoretical value. The agreement of both confirms the loss of contrast for long interferometer times 
2T — ST. For smaller SNRs, the minimum detectable contrast increases, as indicated by the dashed line. 
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Asymmetric Mach-Zehnder interferometer 

We now turn to the theoretical description of our asymmetric Mach-Zehnder interferometer 
(AMZI). Here we only outline the essential ingredients of our formalism and present the main 
results relevant for the experiments described in our article. A more detailed treatment will be 
subject of a future publication. 

Building blocks 

We describe the BEC in our trap by a time- and position-dependent macroscopic wave function 
ip = ip(r,t) determined by the non-linear Gross-Pitaevskii equation (GPE). In all experiments 
reported in our article the trap is well approximated by a parabolic potential V which is also 
time-dependent since the BEC is released by turning off the trap. Another reason for considering 
a time-dependent V emerges from the need to describe the technique of delta-kick cooling in our 
experiments. In these cases, the GPE can be well approximated by an almost analytical solution 
ip = e** ^ TF \r, t) consisting of the product of a phase factor and the Thomas- Fermi (TF) ground 
state wave function of the trap with a time-dependent rescaling of the spatial coordinates^' 32 ~ 3 ^. 
The so-defined time-dependent wave function ip( TF \r, t) is real and normalized at all times. More- 
over, the phase $ of ip does not only depend on the space-time variables r and t but also on the 
center-of-mass (COM) motion of the BEC expressed by the time-dependent coordinate 1Z = TZ-(t) 
and momentum P = P(i). In contrast, the TF-wave function ^;( TF ) is only a function of the 
difference r — TZ.it) and t. 

The interaction of the BEC with a Bragg pulse at a time t p can change the instantaneous 
momentum P(i p ) of the COM motion by ±hk, where k is the effective wave vector of the Bragg 
pulse, but leaves the position 7Z.(t p ) untouched. Moreover, the wave imprints a phase on the BEC 
given by the sum of k ■ TZ-(t p ) and the phase <t>(t p ) of the Bragg pulse. The wave function altered 
in this way is also multiplied by the imaginary factor (— i) resulting from the time-dependent 
Schrodinger equation. 

When the parameters of the Bragg pulse are chosen to act as a 50:50 beam splitter, the BEC 
partly experiences a recoil and partly continues undisturbed with equal amplitudes. As a result, 
the COM motion is in a superposition of the original momentum P and "P + hk, or P and "P — hk 
with equal weights 1/V2. When the pulse acts as a mirror, it changes the momentum by +hk or 
—hk with unit probability. 

Wave functions at the exit ports 



The description of the BEC with the help of the GPE and the effect of the Bragg pulses outlined 
above constitute the main building blocks of our formalism to analyze the AMZI shown in Fig. A5 
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Our goal is to obtain the wave functions and of the BEC at the two exit ports I and II of 
the interferometer, which differ by the COM momentum of the BEC. 



y 
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Figure A5. Diagram of the COM trajectories for our AMZI with labels for the two branches (a and b), the 
two exit ports (I and II) and the equal-momentum trajectories at each exit port ({1,2} and {3,4}). The 
horizontal axis corresponds to time and the vertical one to position along the y-direction. The dashed lines 
depict the Bragg laser pulses whereas the dotted line represents the detection pulse for absorption imaging. 

Indeed, in port I it is the initial momentum "P since both paths 1 and 2 have experienced two 
opposite momentum kicks. However, in port II it is "P + hk since both paths have received an odd 
number of kicks which is different in paths 3 and 4. In path 4 we have three exchanges of momenta, 
that is hk — hk + hk = hk, whereas in path 3 we have only a single one, that is +hk. 

In order to find the wave functions and ^>u at the two exit ports, we now solve the GPE 
including the momentum kicks imparted by the Bragg pulses starting with the initial TF-wave 
function of the ground state of the trap. We switch off the trap at t = and evolve the BEC 
according to the GPE. This dynamics gives rise to a phase factor <3?o and a spreading of the TF- 
wave function. The Bragg pulse at t\ = Tq serving as a beam splitter either changes the COM 
motion due to the momentum transfer hk and imprints phases on the BEC in addition to the $o 
accumulated before the pulses, or leaves it untouched. After the pulse the two BECs evolve for 
the time T according to the GPE in the absence of a potential, which provides us with additional 
dynamical phases. 

The second Bragg pulse at t2 = To + T acts as a mirror and imprints again a phase but also 
transfers momentum with unit probability. Hence, no splitting of the trajectories occurs. The 
ensuing dynamics of the BECs according to the GPE adds new dynamical phases, and the two 
BECs on the paths a and b approach each other till at = Tq + 2T — 5T the third Bragg 
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pulse implementing another beam splitter mixes them. Again the phases of the pulse and the 
instantaneous positions of the BEC are mapped onto the wave functions. 

At this moment it becomes important that our interferometer is asymmetric since the time T 
between the first and the second Bragg pulse is different from the one between the second and the 



third by an amount ST, as indicated in Fig. A5 As a consequence, at each exit port the interfering 
BECs are slightly shifted with respect to each other by an amount d = (hk/m)5T where m is the 
mass of the atom. 

When this displacement d is small compared to the extension of the BEC, we can neglect it 
in the arguments of ip( TF \ which acts as envelope of the wave functions contributing to each exit 
port. However, it is of crucial importance for the phases accumulated by the BEC on the two paths 
through the interferometer, since they depend in a sensitive way on the COM motion. 

An execution of this procedure yields for t 3 < t the expressions 



e i«i^(TF) ( r _7?, 1 (t) ; t) < jl +eX p 



and 



*H = \ e iai1 ^ (TF) (r - n 3 (t);t) jl - exp 
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-fci(t)) 



(6) 



(7) 



for the wave functions \&i and \&n in the two exit ports. Here a\, an and a contain the phases 
imprinted by the Bragg pulses as well as the dynamical phases accumulated by the BEC along the 
two paths and determined by the GPE. 

The two wave functions are centered around the positions 72-i(t) and Ti- 3 {t) of the COM trajec- 
tory at t and the spatially dependent phase, leads to oscillations with wavelength 

2nht 



A = 



m\d\ 



(8) 



BEC densities: Near-field versus jar-field 



From Fig. |A5| we note that at the final beam splitter, that is at £3, we have the equality 
72-1 (£3) = 7^-3(^3)) and due to the momentum difference of hk between TZ 3 and TZ± we find the 
relation 



hk 

K 3 (t 3 + T) =ni{h) + — t 

m 



(9) 



for later times. 

When (h\k\/m)T is larger than the extension of the BEC given by ip( TF \ we can separate the 
two exits and the BEC densities W\ and W\\ following from Q and ([7]) read 



Wi = (V {TF) (r--fci(i);i)) cos : 
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(10) 
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and 



W n = (V TF) (r-^ 3 W;t)) 2 sm 2 
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Three features stand out most clearly in these expressions valid in the far-field limit: (i) The 
envelopes of W\ and Wu are given by the time-dependent TF-profiles (■0*- TF ^) 2 centered at the 
COM coordinates 72-i(i) and "72-3 (t) of the BEC at t. (ii) The arguments of the trigonometric 
functions in W\ and Wu, which are also measured with respect to 7Z-i(i) and 72-3 (i), have identical 
wavelengths given by ^ and depend on the total time from the release of the BEC from the trap till 
the measurement, denoted by T ex in the paper, and (Hi) the two distributions are complementary, 
that is Wj follows a cosine- and Wu a sine-function. 

It is also interesting to note that we would have obtained the same expression for W\ if we had 
considered the free propagation of a superposition of two identical wave packets initially separated 
by d and with momentum "P. Indeed, the fringe spacing A in both situations is identical. In the 
corresponding analogy giving rise to Wu the wave packets have the momentum "P + hk. It is in 
this sense that the AMZI is analogous to Young's double-slit configuration. 

In the opposite limit, when {h\k\/rn)T is smaller than the extension of the BEC, the two wave 
functions and \Pn create an interference term / = ^J^ii + c.c. in the probability 

W = |^i + ^ii| 2 = Wi + Wn + I . (12) 

From ([6]) and ([7]) we note that / contains the phase difference a\ — an, which apart from many other 
phases involves a term k ■ r due to the different COM momenta of the paths 3 and 1 after the last 
beam splitter and gives rise to oscillations of the term / with the optical wavelength. Our detectors 
cannot resolve such small structures and average instead over them, so that the contributions of 
the interference term I vanishes. 

When we now arrange the detection time r such that 

±JE - 1 (i 3 ) 

we can make a cosine- out of the sine- function in Wu, which amounts to a shift by A/2 so that the 
maxima of W\ and Wu coincide. 

Moreover, we recall that we are in the near-field regime where we cannot distinguish the two 
exit ports, that is we can replace the COM coordinate 72- 3 by IZi in the prefactor ?/;( TF ) of Wu, 
which finally yields the expression 



W = 2 (V TF) (r -K^t);^) 2 cos 2 
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(14) 



for the BEC density. 

Hence, by working in the near field region, we have obtained a factor of two in the intensity of 
W which allowed us to obtain a higher signal and extend our measurements to longer times. In 
Fig. 3 such measurements are indicated by triangles. 



